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ABSTRACT 

The problem of construction of a general ihomogeneous solution of /^-dimensional Einstein 
equations with any matter sources satisfying the inequality e > p in the vicinity of a cos- 
mological singularity is considered. It is shown that near the singularity a local behavior of 
metric functions is described by a billiard on a space of a constant negative curvature. The 
billiard is shown to have a finite volume and consequently to be a mixing one. Dynamics 
of inhomogeneities of metric is studied and it is shown that its statistical properties admit 
a complete description. An invariant measure describing statistics of inhomogeneities is ob- 
tained and a role of a minimally-coupled scalar field in dynamics of the inhomogeneities is 
also considered. 
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1 Introduction 



As is well known a number of unified theories predict that dimension of the Universe exceeds 
that of we normally experience at macroscopic level [1]. It is assumed that presently addi- 
tional dimensions are hidden, for they are compactified to the Planckian size, and they do 
not display themselves in macroscopic and even in microscopic processes. However, the situ- 
ation must be changed as we come back with time to the very beginning of the evolution of 
our Universe. Standard cosmological models predict the existence of a singular point at the 
very beginning and, therefore, the universe size could approach to the Planckian scale. Thus, 
in the early universe the additional dimensions, if exist, must not be different from ordinary 
dimensions and should be taken into account. Moreover, one could expect that the existence 
of additional dimensions may drastically change properties of the singularity and even re- 
move it. The main aim of this paper is to construct a general solution of multidimensional 
Einstein equations near a singularity and to investigate properties of inhomogeneities. 

The way to construct a general solution with singularity was indicated first by Belinsky 
et al.in Ref. [2] for D — A, where D is the dimension of a spacetime. Dynamics of metric at a 
particular point of space was shown to resemble the behaviour of the well studied "mixmas- 
ter" (or of the type- IX) homogeneous model and the last one has a complex stochastic nature 
[3,4]. Subsequent utilizing of that construction has been done in Ref. [5] where the so-called 
scalar-vector-tensor theory (or the case D — 5) was considered and the main feature of the 
mixmaster model, i.e. the complex oscillatory regime was shown to be also present in the 
5-dimensionaI case. 

An investigation of inhomogeneities of metric based on the general solutions has been 
considered first in Ref. [6] . The case of the scalar-tensor theory (or D = 4+ scalar fields) 
was considered and it turned out that the oscillatory regime leads to the fractioning of the 
coordinate scale A of the inhomogeneities of Kasner exponents (A ^ Ao2~^, where N is the 
number of elapsed Kasner epochs and Aq is the initial scale of inhomogeneities). However, 
the methods by means of which the properties and statistics of the inhomogeneities were 
investigated turned out to be unapplicable for general case (i.e. for the absence of scalar 
fields as well as for the expanding universe). This problem has been solved recently in Ref. 
[7]. In this paper we generahze the results obtained in Ref. [7] to the case of arbitrary 
number of dimensions D. 

As it was mentioned above the main features of the dynamics of an inhomogeneous 
gravitational field nearby the singularity in 4-dimensional case may be summarized as follows: 

1. Locally dynamics of metric functions resembles the behaviour of the most general ho- 
mogeneous "mixmaster" model [2], which has stochastic behaviour [3,4]. Just the stochastic 
behaviour leads to a monotonic decrease of the coordinate scale of the metric inhomogeneities 
[6,7]. 

2. In the vicinity of a singularity a scalar field is the only kind of matter effecting the 
dynamics of metric [5]. 

These facts may be simply understood under the following qualitative estimates (that 
is confirmed by subsequent consideration). As is well known in cosmology the horizon size 
Ih is a natural scale measuring a distance from the singularity. Therefore, inhomogeneities 
may be divided into the large-scale (/j ^ //j) and small-scale (/j -C Ih) ones. The horison 
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size varies with time as Ih t (where t is the time in synchronous reference system) whereas 
the characteristic spatial dimension of the inhomogeneity may be estimated as ~ (as 
t ^ 0). In a hnear theory for an isotropic background the exponent a may be expressed 
via the state equation of matter as a; = 3(|^ and what is important a < 1. Thus, it is 
clear that an arbitrary inhomogeneous field becomes large-scale in the sufficient closeness to 
the singularity. Since the inhomogeneities are large-scale there are no effects connected with 
propagating of gravitational waves etc, and this would mean that inhomogeneities become 
passive. Consequently, dynamics of the field may be approximately described by the most 
general homogeneous model depending parametrically upon the spatial coordinates. Note, 
however, that the homogeneous model would appear to be in a general non-diagonal form. 

The second fact may be understood in the same way. As it was shown in Ref. [8] the 
gravitational part of the Einstein equations at the singular point varies with time, in the 
leading order, as Rp ~ whereas the matter has the order TJ' ~ t"^*^, where k depends 
upon the state equation as /c = Thus, one can see that for the equation of state satisfying 
the inequality p < e we have k < 1 and only for the limited case p = e the both sides turn 
out to be of the same order. We note that in the vicinity of a singularity scalar fields give 
just this equation of state. 

As it is well known (see for example Ref. [5,9]) additional dimensions may be treated in 
ordinary gravity as a set of nonminimally coupled scalar and vector fields. Therefore, one 
could expect that the main contribution in dynamics in the vicinity of a singularity would 
be given by those dynamical functions which are connected with scalar fields whereas other 
functions would play a passive role. 

Thus, one could expect that in multidimensional cosmology local behaviour of the metric 
functions (at a particular point of space) will be described by a most general homogeneous 
model. Here, it is necessary to recall the important property of the mixmaster universe that 
is the stochastic behaviour. The problem of stochasticity of homogeneous multidimensional 
cosmological models has been investigated in a number of papers [10]. In particular, the 
conclusion was made that chaos is absent in spaces whose dimension is large enough. This 
negative result is, apparently, connected with the fact that in the vicinity of a singularity 
the homogeneous model would appear to be in a general nondiagonal form, instead of simple 
diagonal models considered in Refs.[10]. Besides, recently in Ref.[ll] it was shown that a 
wide class of multidimensional models has stochastic behaviour. In this paper we shall show 
that stochastic behaviour is the general property of multidimensional cosmology. 

2 Generalized Kasner Solution, Generalized Kasner Vari- 
ables 

We consider the theory in canonical formulation. Basic variables are the Riemann metric 
components Qa/s with signature (+, — , — ) and a scalar field specified on the n-manifold 
S, and its conjugate momentum 11"^ = ^Jg{K"f^ — g"'^K) and U^, where a = 1, ...,n and 
K"'^ is the extrinsic curvature of S. For the sake of simplicity we shall consider S to be 
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compact i.e. dS — 0. The action has in Planck units the following form 



+ - TVi/O - No,H'')(F'xdt, (2.1) 



where 



here 



^° = ^ { n^n^ - ^^{Kr + + 9{w{<p) -r)}, (2.2) 
i/° = -2n|^/ + ^-^9^0n^, (2.3) 

= I { 9''^do,(t>M + V{(l>) } . (2.4) 
A generalized Kasner solution is realized under the following assumption 

~ (n^, n^) » y = g{w - r), (2.5) 

where ^/gT denotes the first three terms in (2.2). Then, using (2.1) one can find the following 
solution of the multidimensional Einstein equations 

n— 1 

ds^ ^dt^-J2 i'^ll > l^dx'^dx'^ (2.6) 

a=0 

where Sa are fimctions of space coordinates. Kasner exponents Sa satisfy the identities 

= Y,s\ + = I, and run the domain — < Sa < 1 (here = 2 (n^V ^^' 
as it was shown in Ref. [2,5] the generalized Kasner solution takes a substantial portion of 
the evolution of metric it is convenient to introduce a Kasner-like parametrization of the 
dynamical variables [7]. We consider the following representation for metric components 
and their conjugate momenta 

9ap^J2^w{q^}lllh (2-7) 

a 

^p = EPaL:ih (2.8) 

a 

here L^l^ — 6^ {a,h — 0, {n — 1)), and the vectors contain only n{n — 1) arbitrary 
functions of spatial coordinates. Further parametrization may be taken in the following 
form 

II = U^Sl VI e SOin), S^^ = 5: + Rl (2.9) 

where i?^ denotes a triangle matrix (i?^ = as a < a). Substituting (2.7) - (2.9) into (2.1) 
one gets the following expression for the action functional 

/ = j^{Pa^ + + - NH' - N^H^)d^xdt, (2.10) 

here — 2J2bPbL^U^ and the Hamiltonian constraint takes the form 

= ^ { EpI - ^^iJ:Par + m + V}. (2.11) 



In the case of n = 3 the functions are connected purely with transformations of a 
coordinate system and may be removed by solving momentum constraints if" = 0. In the 
multidimensional case the functions i?^ contain "^""^-^ dynamical functions as well. Now it 
is easy to see that the choice of Kasner-like parametrization simplifies the procedure of the 
constructing of the generalized Kasner solution. Indeed, if we now neglect the potential term 
in (2.10) and put N"' = we find that Hamiltonian does not depend on the scale functions 
and other dynamical variables contained in Kasner vectors introduced by expressions (2.7) 
(2.8). 



3 The asymptotic model in the vicinity of a cosmolog- 
ical singularity 

As it is well known, [2], [5], the Kasner regime (2.6) turns out to be unstable in a general 
case. This happens due to the violation of the condition (2.5) because the potential V 
contains increasing terms which lead to replacement of Kasner regimes. To find out the law 
of replacement it is more convenient to use an asymptotic expression for the potential [7], 
[11]. For this aim we put the potential in the following form 

k 

V=Y.^Ag''\ (3.1) 

A=l 

here Aa is a set of functions of all dynamical variables and of their derivatives and Ua are 
linear functions of the anisotropy parameters Qa = (ma = ua{Q))- Assuming the 
finitencss of the functions A and considering the limit g ^ we find that the potential V 
may be modeled by potential walls 

9"-«.MQ)l = {ot~'T>°o 

Thus, putting N°' = we can remove the passive dynamical function T^, from the action 
(2.10) and get the reduced dynamical system 

I- I {Paf^ + n^f -A{Ep'-^(Ep)' + in^ + c/(g)}}d"xdi, (3.3) 

here A is expressed via the lapse function as A = In harmonic variables the action (3.3) 
takes the form formally coincided with the action for a relativistic particle 

7 = /■ { P,^ -X{Pf + U- Pi) } (Txdt, (3.4) 
J s 



here r = 0, i = l,...,n, q" = Ap^ + {j = l,...,n - 1), ^" = V 2^*^ the 
constant matrix Aj obeys the following conditions 

Y: a • =0,Y: A'^Al = n{n - 1)6^, (3.5) 
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and can be expressed in the following form 



n(n — 1) 



where I ]:^.^^ . 
^ i 0,j < a 

Since the timelike variable varies during the evolution as ^ \ng the positions of 
potential walls turn out to be moving. It is more convenient to fix the positions of walls. 
This may be done by using the so-called Misner-Chitre like variables [11] {y = y^) 



Using these variables one can find the following expressions for the anisotropy parameters 

Q.fe) = i{i + ^}. (3.7) 

which are now independent of timelike variable r. From (3.7) one can find the range of the 
anisotropy functions — < Qa < 1- 

Choosing as a time variable the quantity r (i.e. in the gauge N — ^^^^^^y/^exp(— 2r) / P°) 
we put the action (3.4) into the ADM form 

/ = { Pfj+ P--§^z- - P%P, y) } (TxdT, (3.8) 

where the quantity 

P^P^y) = {e\y,P)+V[y] + {P^^fe-^^f/\ (3.9) 
plays the role of the ADM Hamiltonian density and 

The part of the configuration space connected with the variables y is a realization of the 
(n — 1) -dimensional Lobachevsky space [12] and the potential V cuts a part of it. Thus, 
locally (at a particular point of S) the action (3.9) describes a billiard on the Lobachevsky 
space. The positions of walls which form the boundary of the billiard are determined, due 
to (3.1) by the inequalities 

<Jabc = l + Qa-Qb-Qc>^.a^h^c (3.11) 

and the total number of walls is "("--^K^-^) ^ Usj^g matrix (3.5) one can find that the 
walls are formed by spheres determined by the equations 

<^abc = { {y + B^bc? + 1 - > B^^c = :^^{A'^ A% (3.12) 



here for arbitrary a, 6, c we have B'^ — 1 + In a general case n points of the biUiard 

— * — * 

having the coordinates Pa — he on the absolute (at infinity of the Lobachevsky space). 

Nevetheless, one can show that the volume of the billiard is finite. We give two simplest 
examples for illustration of the billiards on fig.l. The case n = 3 on fig. la coincides with the 
well-known "mixmaster" model and on fig. lb we illustrate the case of n = 4 considered in 
Ref.[5]. 



4 Dynamics of inhomogeneities 

The system (3.8) has the form of the direct product of "homogeneous" local systems. Each 
local system in (3.8) has two variables e and as integrals of motion. The solution of this 
local system for remaining functions represents a geodesic fiow on a manifold with negative 
curvature. As it is well known the geodesic fiow on a manifold with negative curvature is 
characterized by exponential instability [12]. This means that during the motion along a 
geodesic the normal deviations grow no slower than the exponential of the traversed path 
^ ~ Coc*) , where the traversed path is determined by the expression 

This instability leads to the stochastic nature of the corresponding geodesic fiow. The system 
possesses the mixing property [13] and an invariant measure induced by the Liouviullc one 

dn{y, P) = const5(E - e)d''-^yd''-^P, (4.2) 

where is a constant. Integrating this expression over e we find 



dl^iV: s) = const— (4.3) 



where s——, Isl = 1 



e 



Since the inhomogeneous system (3.8) is the direct product of " homogeneous" systems one 
can simply describe its behaviour as in ref [7]. In particular, the scale of the inhomogeneity 

decreases as 

A.-(|^)-'-A°exp(-5) (4.4) 

and after sufficiently large time (s(r) — >• oo) the dynamical functions y{x), P{x) become a 
random functions of the spatial coordinates. Their statistics is described by the invariant 
distribution (4.3) and asymptotic expressions for averages and correlating functions have the 
form 

< y{x) >=< P{x) >= 0, < yk{x),yi{x') >=< yk, yi > 5{x, x'), (4.5) 

for \x — x'\ 3> A^exp(— s). 

Here it is necessary to point out a role of the scalar field in dynamics and statistical 
properties of inhomogeneities. As may be easily seen from (4.1) in the absence of a scalar 
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field (i.e. P" — 0) the transversed path coincides with the duration of motion (we have 
s — At = r — To instead of (4.1)). Thus, the effect of scalar fields is displayed in the 
replacement of the dependence for transversed path of time variable and, therefore, in the 
replacement of the rate of increasing of the inhomogcncities. This replacement does not 
change qualitatively the evolution of tlic universe in the case of cosmological expansion. 
But in the case of the contracting universe the situation changes drastically. Indeed, in 
the limit r —oo from (4.1) we find that the transversed path s takes a limited value Sq 
and therefore the increasing of inhomogeneities turns out to be finite. One of consequences 
of such behaviour is the fact that at the singularity the functions y and P take constant 
values. In other words in the presence of scalar fields a cosmological collapse ends with a 
stable Kasner-hke regime (2.6). This fact may be seen in the other way. Indeed, in the limit 
T — oo the scalar field gives the leading contribution in ADM Hamiltonian (3.9) and P° 
does not depend on gravitational variables at all. 

The finitcncss of the transversed path s(r) leads, generally speaking, to the destruction 
of the mixing properties [13], since for establishment of the invariant measure it is necessary 
to satisfy the condition sq ^ oo. Evidently, this condition requires the smallness of the 
energy density for scalar field as compared with the ADM energy of gravitational field (the 
last term in (3.9) in comparison with the first ones). Indeed, in this case sq is determined 
by the expression sq = — In , which follows from (4.1), and as P" one get sq — >■ oo 
(i.e. s can have arbitrary large values). 

Thus, in the case of cosmological contraction one may speak of the mixing and, therefore, 
of establishment of the invariant statistical distribution just only for those spatial domains 
which have suficiently small energy density of the scalar field. 

5 Estimates and concluding remarks 

In this manner the large-scale structure of the space in the vicinity of singularity acquires a 
quasi-isotropic nature. A distribution of inhomogeneities is determined by the set of functions 
of spatial coordinates e(x), n^(x) and which conserve during the evolution a primordial 
degree of inhomogeneity of the space. The scale of inhomogeneity of other functions grows 
as A ~ Aoe^**^^). In this section we give some estimates clarifying the behaviour of the 
inhomogeneities. For simplicity we consider the case when the scalar field is absent. 

To find the estimate for growth of the inhomogeneity in a synchronous time t {dt — Ndr) 
we put y — Then for variation of the variable r one may find the following estimate 

^ cxp(— ^e"'^) ~ P°t, (here the point t = corresponds to the singularity). According 
to (4.4) the dependence of the coordinate scale of inhomogeneity upon the time t takes the 
form 

Aohi(l/^o)/ln(l/<7) 

in the case of contracting {g — > 0) and 

Af^ Aoln(l/^)/ln(l/^o) 
in the case of the expanding universe. 



7 



A rapid generation of the more and more small scales leads to the formation of spatial 
chaos in metric functions and so the large-scale structure acquires a quasi-isotropic nature. 
Speeds of the scale growing (Hubble constants) for different directions turn out to be equal 
after averaging over a spatial domains having the size a; Aq. Indeed, using (3.7) one may 
find the expressions for averages < Qa >= l/n. 

Besides, it is necessary to mention one more characteristic feature of the oscillatory 
regime in the inhomogeneous case. This is the formation of a cellular structure in the scale 
functions Qa during the evolution which demonstrate explicitely the stochastic process of 
development of inhomogeneities. Indeed, let us consider some region of coordinate space 
AV. Two functions y{x) define the map of that region on some square T, ^ K (see fig.lc). 
During the evolution the size of the square E grows ~ e**^"^) and E covers the domain of the 
billiard K many times. Each covering determines its own preimage in AV. In this manner 
the initial coordinate volume is splitted up in "cells" AV = [j^AVi. In the every cell the 
vector y{x) takes almost all admissible values y E K and that of the functions Qa (for n = 3 
Qa e [0, !])• Such a structure turns out to be depending of time and the number of cells 
increases as N Noe^^'^\ However, the situation will be changed if we consider a contracting 
space filled with a scalar field. Then the evolution of this structure in the hmit g ^ ends, 
because the functions Qa become independent of time, and on the final stage of the collapse 
one would have a real cellular structure [6]. 

In spite of the isotropic nature of the spatial distribution of the field the large local 
anisotropy displays itself in the anomalous dependence of spatial lengths upon time vari- 
able for vectors and curves . Indeed, in the case oi D — 3 a moment of scale function 
< (yf^^*3" > (where M > 0) is decreased in the asymptotic o as the Laplace integral 
!o 9^^°^ p{Qa)dQa, where p{Qa) is a distribution which follows from (4.3) and has the form 



As Q <^ 1 one has p{Qa) ~ f (Qa) ^^'^ and, thus, in the limit g' — > we get the estimate 



When D > 3 the similar analysis may be done. 

In conclusion we briefiy repeat the main results. The general ihomogeneous solution of 
^^-dimensional Einstein equations with any matter sources satisfying the inequality e > p 
near the cosmological singularity is constructed. It is shown that near the singularity a local 
behavior of metric functions ( at a particular point of the coordinate space) is described by a 
billiard on the (D — l)-dimensional Lobachevsky space. In contrast to diagonal homogeneous 
models [10] in the inhomogeneous case the billiard is shown to have always a finite volume 
and consequently to be a mixing one. The rate of growth of inhomogeneities of metric is 
obtained. Statistical properties of inhomogeneities are described by the invariant measure. 
It is shown that a minimally-coupled scalar field leads, in general, to the distruction of 
stochastic properties of the inhomogeneous model. 
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<g^«« >^ (Mln(l/y))-V2. 



(5.2) 
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